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Abstract 

Normality in connection with 75-hermiticity determines the basic chiral proper- 
ties and rules. The Ginsparg- Wilson (GW) relation is one of the allowed constraints 
on the spectrum. Interrelations between features of the spectrum, the sum rule for 
chiral differences of real modes and contributions to the Ward identity are pointed 
out. The alternative chiral transformation of Liischer gives the same Ward identity 
as the usual one, in the global and in the local case. Imposing normality on a gen- 
eral function of the hermitean Wilson-Dirac (HWD) operator, inevitably leads at 
the same time to the Neuberger operator and to the GW relation. In this context 
also the case with zero eigenvalues of the HWD operator is handled. The eigenvalue 
flows of the HWD operator obey a differential equation the characteristic features 
of the solutions of which are discussed. 
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1. Introduction and overview 



Recently considerable progress in the description of chiral fermions has been initiated by 
works of Neuberger [Q] from the point of view of the overlap formalism and of Hasenfratz 
et al. in the context of fixed point actions. This has also revived considerations of the 
Ginsparg- Wilson (GW) relation |^ which in both cases turned out to be satisfied. On 
the basis of the GW relation Liischer ||^ proposed an alternative chiral transformation 
providing chiral symmetry at the classical level. The finite form of this transformation 
has been given by Chiu 0. 

Neuberger in particular, was able to derive an explicit form of the massless Dirac 
operator based on the hermitean Wilson-Dirac (HWD) operator, which also plays a major 
role in the overlap formalism These developments have then given rise to several 
numerical studies of eigenvalue flows of the HWD operator with the mass parameter 0, 
relying on the fact that this operator has well defined spectral properties. 

Despite the many publications which followed the mentioned works, there are clearly 
still many questions open. In the present paper we address a number of them which are 
related to a more precise understanding and to basic properties of these new developments. 

We start from the observation that the Dirac operator must be normal in order that 
reliable conclusions become possible. This follows from two theorems of the spectral theory 
of operators in unitary spaces. The first one of them says that normality is necessary and 
sufficient in order that the eigenvectors form a complete sytem. Apart from the fact that 
otherwise very little is known on spectral properties, this implies that without normality 
there are necessarily defects which can hardly be tolerated in making predictions. 

The second one of these theorems more specifically concerns chiral properties. It states 
that normality guarantees that an eigenvector of the operator is at same time an eigen- 
vector of the adjoint operator. This together with 75-hermiticity (i.e. hermiticity of the 
operator multiplied by 75) is exactly what provides the basis for chiral behaviors. 

We point out that, given normality and 75-hermiticity, one already obtains the basic 
rules and properties. Apart from the general structure of the eigenvector system, the 
fundamental sum rule for chiral differences of real modes emerges. The relations between 
the operator and such modes, and in particular those to the index, as needed in various 
contexts, immediately follow. 

A further consequence of normality is that the operator decomposes into commuting 
operators related to the real part and to the imaginary part of its eigenvalues. This 
allows to study possible constraints on the spectrum. In particular, restricting it to a 
one-dimensional set, we find that, in addition to zero, the curve must meet the real axis 
at least at one further point, in order that the sum rule mentioned above allows a nonzero 
index. From this point of view the general nature of the GW relation becomes clear. It 
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is just one of such constraints which satisfies the requirement of a further real eigenvalue 
in a minimal way. 

The symmetry at the classical level provided by the alternative chiral transformation 
1^, makes things similar to what one is accustomed to in continuum theory, however, 
at the price of complications due to the action dependence of this transformation. At the 
quantum level, where there is anyway no such symmetry, the question arises, what the 
precise difference to the usual chiral transformation is. 

In order to have a basis for a general comparison of transformations, we derive the Ward 
identities in an appropriately general way. It is seen that by the normality of the operator 
the global chiral transformation leads just to the sum rule for chiral differences of real 
modes. The correspondence of the contributions to terms familiar in continuum theory is 
pointed out. The results of the local chiral transformation are similarly identified. 

It is shown that the alternative chiral transformation leads exactly to the same result as 
the usual chiral one, even without assuming the GW relation. Imposing the GW relation 
has the only effect to specialize the results to that case with only two real eigenvalues. We 
also give the local version of the alternative transformation. The resulting Ward identity 
is again seen to agree with that of the usual local chiral one. 

The operator of Neuberger is the only explicit form of a massless Dirac operator on 
the lattice presently known. In order to get more insight into the possibilities of the 
construction of such operators it appears desirable to have also a derivation of it which 
does not rely on the overlap formalism. Further, there is the somewhat unsatisfactory 
point that zero eigenvalues of the HWD operator so far had to be excluded. 

In our derivation of the indicated operator the requirement of normality is central. To 
avoid doublers and at the same time to deal with well defined mathematical properties, 
at present the only possibility is to start from the HWD operator. Therefore we consider 
a general function of this operator and impose the necessary conditions on it. Doing this 
it turns out that the requirement of normality is an extremely strong one, leading at the 
same time to the Neuberger form of the Dirac operator and to the GW constraint on the 
spectrum. In addition the inclusion of zeros of the HWD operator gets nonstraightforward. 
Nevertheless, also for this problem a way out is found. 

For the explicit Dirac operator the relations of its eigenvectors to that of the HWD 
operator become transparent. The apparent importance of the HWD operator suggests 
to study its eigenvalue fiows with the mass parameter also analytically. We show that they 
satisfy a differential equation and give a complete overview of the characteristic properties 
of its solutions. 

Section 2 is devoted to the basic chiral properties. In Section 3 possible constraints on 
the spectrum are discussed. Section 4 contains the general derivation of Ward identities. 
In Section 5 the results for particular transformations are analyzed. Section 6 gives the 
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systematic construction of a normal operator. In Section 7 the spectral flows of the HWD 
operator are investigated. 



2. Basic chiral operator properties 

We require D to be normal 

[A ^1=0 (2.1) 

and 75-hermitean 

^^ = 75^75 • (2.2) 
Because of ( p.l|) the solutions fk of the eigenequation 

Dfk = Xkfk (2.3) 

form a complete orthonormal set, on the basis of which general conclusions become pos- 
sible. By (|2.1|) simultaneous eigenvectors of D and exist. In the present context this 



has the important consequence that one also has 

D^fk = Kfk (2.4) 

which together with ( f2.2D leads to 

Di,fk = Xll,fk . (2.5) 



The comparison of ( |2.3| ) multiplied by 75 with (|2.5|) then gives 



[j,,D]fk = if Aureal, (2.6) 

which tells that in the subspace of real eigenvalues of D one can introduce simultaneous 
eigenvectors of D and of 75, i.e. ones with chirality. 
Multiplying (|2.3|) from the left by 7/75 and its adjoint fiD^ = flX*i from the right by 



75/fc one obtains the relation 

//75/fe = for , (2.7) 

which actually reflects the orthogonality of eigenvectors related to different eigenvalues of 
D. In view of ( |2.6|) and of the comparison of ( p.5|) with (|2.3| ), respectively, it is convenient 
to introduce in more detail 

/f for ImAfc = with 75/f = c^/f 
fk={ fi:^ for ImAfc>0 with 75/f = /f (2-8) 
/f for ImAfc<0 with 75/f = 
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where for the chirahty Cr possible values are +1 and —1. Obviously eigenvectors related 
to complex eigenvalues always come in pairs, while those related to real ones need not to 
do so. Of course, how many of each type occur depends on the particular D, however, by 
and ( |2.2| ), the structure of the eigenvector system is in any case the one described. 
For the numbers of modes related to a real eigenvalue A of D 



1 ± Cfc 



(2.9) 



A f. — A real 



from (|2.7|) , (p.8|) and Tr(75) = one obtains the sum rule for the chiral differences of real 
modes 

E (iV4A)-iV„(A))=0 . (2.10) 



A real 



It in particular implies that the difference for eigenvalue zero A^_(0) — A^+(0), the index of 
D, can only be nonvanishing if a corresponding difference from nonzero eigenvalues exists. 

From ( p.7|) , (p.8|) and ( p.9|) one also readily gets the useful relations 



£Tr((D + e)-^-f5) N+{0) - N_{0) for e^O 



Tt{{D 



I5D) 



E 

\=iO real 



iV+(A) -iV_(A) 



for 



Tr(75l^)= A(iV+(A)-iV_(A)) 

At^O real 



(2.11) 
(2.12) 

(2.13) 



between D and A^+(A) — A^_(A). It will be discussed later that ( ^.ll] ) and ( |2.12| ) are related 
to the index term and the topological charge (FF) term, respectively, of continuum theory. 
They obviously add up to the sum rule (|2.10|) . If only one nonzero value for real eigenvalues 
occurs, one can also use (|2.13|) instead of (|2.12|) to relate D to the corresponding chiral 
difference. 



3. Constraints on location of spectrum 



In order to study possible constraints on the spectrum of D we use the decomposition 



D = u + iv with u = u^ = ^{D + D^) 



V = v 



1 



{D-D^ 



for which (|2.1|) , the normality of D, implies 

[m, = 



(3.1) 



(3.2) 



By ( p.2|) , the 75-hermiticity of D, it follows that ( |3.1| ) is at the same time the decom- 
position into the parts commuting and anticommuting with 75 and that [75, m] = and 
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{75; v} = hold. According to (|3.2|) u, v and D have simultaneous eigenvectors. The real 
eigenvalues of u and v are simply the real and imaginary parts, respectively, of those of 
D. Therefore we can specify the location of the spectrum of D by constraining u and v. 

A particular simplification arises if the spectrum is located on a one-dimensional set. 
This can be realized by imposing the condition 

J^{u,v) = (3.3) 

with a suitably chosen function J-'{u, v) which, as a function of commuting hermitean op- 
erators, is well defined. In the notation of (p.3|), it satisfies the eigenequation J^{u,v)fk = 
jF(ReAfc, ImAfc)/fc and has the spectral representation J^{u,v) = X^fc -^(ReAfc, ImAfc)/^/^. 
Because we wish to allow for the eigenvalue of Z^, the function considered as a func- 
tion of real numbers should have the property JF(0, 0) = 0. Since we in addition want 
that the index of D can be nonzero, according to ( |2.10| ) also at least one real eigenvalue 
different from zero must be possible. Therefore in addition one has to require JF(/3, 0) = 
for at least one real (3^0. Thus the function must have the particular properties 

J^(0, 0) = and 0) = for some p (3.4) 



in order that the curve specified by ( |3.3| ) meets the real axis at zero and at least at one 
further point. 

A simple possibility which satisfies condition ( p.4| ) is a circle through zero with center 
on the real axis 

J^{u,v) = {u- pf + - p'^ = . (3.5) 
Using ( |3.1D and ( p.2| ) one can write ( |3?5| ) as p{D + D"^) = D^D which by ( pl2| ) is seen to 



be just the GW relation 

{^,,D} = p-'D^,D . (3.6) 



In contrast to the original form [Q, however, (|3.6| ) does not involve a further operator 
sandwiched in addition to 75 into its r.h.s. because this would spoil the normality of D. 
Therefore here only a real constant p~^ remains. Obviously ( p.5| ) meets the requirement 
( ^.4[ ) in a minimal way, admitting only 2p in addition to 0. The sum rule ( ^.l(j| ) then 
simplifies to two terms. This conforms with the observation of Chiu that in the case 
of the GW relation {75, D} = Dj^D the chiral differences obtained at and at 2 add up 
to zero. 

By ( |3.4| ) the choice J-" = u, corresponding to {75, D} = 0, is excluded because there 
is only one real eigenvalue (on the finite lattice with D being bounded a second one 
can also not occur at 00). This choice may, however, be approached in the continuum 
limit. In fact, considering the stereographic projection of the circle ( p.5[ ) on the sphere of 
complex numbers it is seen that for increasing radius it approaches the circle through 00 
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on this sphere which corresponds to the imaginary axis in the plane. With p ~ 1/a for 
decreasing lattice spacing a the envisaged approach indeed occurs. This suggests that in 
the continuum the sum rule for chiral differences could possibly be satisfied by eigenvalues 
at and at oo. Of course, the subtleties of the respective limit remain to be investigated. 

The form of the constraint on the spectrum depends on the particular properties of the 
operator D considered. The Neuberger operator is tied to the GW constraint. For other 
constructions in any case also ( p.4|) is to be required. In addition, also the appropriate 
behavior in the limit should be guaranteed. 



4. General form of the Ward identity 

In dealing with Ward identities it should be remembered that the expectation values 

involve integrals J[dipdip]e~^^ = detM and 

J [d^d^]e-^^^,,^,, . . . = E ^::tM-l, ■ ■ ■ det M (4.2) 

(where ef^" ff = +1, — 1 or if fci . . . kg is an even, odd or no permutation of /i . . . Ig). 
Thus, in order that the expectation values (^?T|) are properly defined, must exist and 
det M be nonzero. Therefore, to be able to proceed in the presence of zero modes of D 
one has to put M = D + e and let e go to zero in the final result. 

Fermionic Ward identities arise from the condition that J[dipdip]e~^^0 must not change 
under the transformation 

ip' = exp{iriT)ilj , ip' = exp^irjT) , (4.3) 

where rj is a parameter. This means that one gets the identity 

A J [d^'d.p']e-'iO'l^^ = (4.4) 

with three contributions, one from the derivative of the integration measure, one from that 
of the action and one from that of O. For the measure contribution within [d'?/''d'?/''] = 
[dV^d?/'] det exp(i?7r) det exp{iriT)^ one obtains 



d (detexp(ir7r) - d {det exp{ir]T) 

drj v=o 0.7] 



= iTrT . (4.5) 

ri=0 ^ ' 
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With respect to the derivative of O we note that for products V oiifj and il) fields, since the 
Grassmann-even com 
place, one can relate 



Grassmann-even combinations V'jg^ and V'j^f^ '^^^ readily shifted to the appropriate 



dP' 



. («) 

The fermionic part of O in general is made up of such products and of sums thereof. More 
specifically it can even be considered to be made up of products of equal numbers of if) 
and if) fields because only such products contribute to the intergrals; for the same reason 
O can also be considered to be Grassmann-even. Thus in any case (|4.6| ) applies to O and 
( [4.4|) becomes 

[d'lpme-^'i -TiiV + V)0 -^iVM + MV)^0 + ^V^ - —V^) =Q . (4.7) 

^ oil) dil> ^ 

This generalizes the relation ||^ which, with suitable choices of O, is used in many appli- 
cations to the case where the integration measure is not invariant. 

In studies of the singlet axial vector current and its relation to anomaly, index and 
topological charge it usually suffices to consider the case O = 1, as is e.g. also done in 
Ref. [||. To keep things general we avoid this here, integrating out the ip and ip fields 
in the second term of ([4.7| ) without specifying O. Then at the same time that term 
gets on equal footing with the first one, as is desirable for a convenient comparison of 
transformations. To integrate out the indicated fields we use the identity 

= l/|d«,)|((AAr'),(e--.V..O) + (Ar'|.)^(e-^*,0)) 

= / |d^.d*|e-a {i,,i„0 + A/^'O + i(^Af-') - (^'S) 

which relies on the fact that j[(\.il)([ip\{d / dipi)G = and J[dil!dip]{d/dil!i)G = for any 
function G. Then (14. 71) becomes 



iW j [dtpdtp]e-^'0 + -J[d'ipdxlj]e-^'[—M-^Rtjj + iljRM-^-^) = (4.9) 

where R = TM - MT and 

W = Ti(-T ~T + M-\TM + MT)) . (4.10) 



To evaluate the terms with derivatives of O in (|4.9|) further we remember that the 
fermionic part of O can be considered to be made up of products of type V = ipj^ipki ■ ■ ■ 
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for which by 



we find 



r - dV r - - dV 

j [di)di)]e-^'—M-^Ri) = + J [dijdij]e-^'ijRM-^-^ 

s E . . . M-\,^_,{M-'RM-'),j^ det M 



(4.11) 



This shows that the terms in (14.91) with derivatives of O cancel and we remain with 



iW J [dipdij]e-^'0 = 



(4.12) 



which inserted into (|4.1| ) gives the Ward identity (WO) = or, if desired, also the one in 
a background gauge field W{0)u = 0. 



5. Results for particular transformations 



For the global chiral transformation, which in terms of ( [4.3[ ) is given by 

r = r = 75 , (5.1) 

the measure contribution — Tr(r + F) vanishes and one obtains 

VF = Tr(M-i{75,M}) (5.2) 



or inserting M = D + e 



W = Tt[{D + £)-H75, D}) + 2eTi[iD + e)-'^. 



(5.3) 



The first term in ( p. 3] ) by ( |2.12| ) is seen to become the sum over 2(A^_|_(A) — A^_(A)) for 
nonvanishing real A and the second one by ( |2.11| ) the difference 2(A^_|_(0) — A^_(0)). By 
they add up to 



PF ^ 2 E {n+W - N4X)) for e 







(5.4) 



A real 



Thus we obviously arrive just at the sum rule for chiral differences of real modes ( p.lO|) 



Analoguous to the features known in continuum theory for quite some time [T^ , the first 
term in (|5.3|) is the the topological charge (FF) term while the second one is the index 
term. The latter by ( p.ll|) is obvious. For the first term the limit has been established long 
ago |]lT| for the Wilson-Dirac operator and recently |jl2[ also for the Neuberger operator 
(for which the use of Tr(75D) by ( |2.12| ) and ( p.l3| ) with only one nonzero A is equivalent 
to using Tr((D + £)-H75,^})). 
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Next we consider the alternative transformation [|, ^ for which in our notation 



r = 75(1 - (2p)-iM) , r = (1 - (2p)-^M)75 . (5.5) 

With ( |5.5| ) one now gets — Tr(r + F) = +p~^TT{'y^M) for the measure contribution and 
Tt{M-\TM + MY)) = Tr(M-H75, M}) - p-^Tii^f^M) for the action contribution. Ob- 
viously the extra term of the latter cancels the measure term so that again the result ( |5.2| ) 
is obtained, and notably, even without assuming the GW relation. 

If with the alternative transformation in addition the GW relation ( p.6|) is imposed, 
inserting M = D + e, for the action contribution one gets Tr((M~^(rM + MT)) = 2e{l + 
{2p)-'^e)TT{{D + ey^'y5) which by ( pTTTD becomes 2(iV+(0) - A^_(0)). For the measure con- 



tribution one has — Tr(r + r) p ^Ti^j^jD) which by ( p.l3| ) equals J2\^o real ^i^+W ~ 
N_{X)). In the GW case with only X = 2p this becomes 2(A^+(2p) - A^_(2p)). Taking 
both contributions[| together it is obvious that one gets the same results as before now 
specialized to the case where only and 2p occur for real eigenvalues. 

The local chiral transformation in the present context is conveniently introduced by 

r = r = 75e(n) with (H^))^,,^, = ^n"Jnn' (5.6) 

for which (^4.10 ) gives 



ly = Tr((M-i{75e(n),M}) . (5.7) 

By inserting the decomposition M = |(M — 75M75) -|- |(M -|- 75M75) (into parts anti- 
commuting and commuting with 75) and also M = D + e into {756(71), M} this becomes 

W = iTr((M-i[e(n), [75, D] ]) + ^Tr((M-i{e(n), {75, D}}) +2eTT{{M-'j,e{n)) . (5.8) 



The first term in ( |5.^ ) is seen to vanish upon summation over n and accordingly corre- 
sponds to the divergence of the singlet axial vector current. Summation over n in the 
rest, responsible for current nonconservation, leads to the results of the global transfor- 
mation. The second term in ( p. 81 ) in the limit gives the FF-density of continuum theory 
]TT| , |T^ . The third term in ( |5.8| ) is the local version of the index contribution. To visualize 
things in terms of current expressions one should remember that the factors in (|5.8|) 
correspond to the integrated out ip and fields. 

We note that the local transformation related to the alternative chiral transformion (|5.5|) 
can also be introduced. It is given by 

r = 75e(n)(l-(2p)-iM) , f = (1 - (2p)-^M)75e(n) . (5.9) 

^In Ref. 1^ the action contribution is missing. The evaluation of the measure contribution there is 
actually a use of the identity = Tr75 = eTr((£) + e)~^75) + Tr{{D + £)^^^c,D) which by inserting the 



GW relation becomes = £(1 + £(2p)-i)Tr((7:> + e)-^-iz) + {2p)-^Tt{-ibD 



10 



The calculation of W with this transformation again leads to (|5.7|) so that it becomes 
obvious that also in the local case nothing new is obtained. 

Of course, also other transformations could straightforwardly be considered along the 
present lines. For the nonsinglet chiral one with flavor operator Ti one has F = F = 'j^Ti in 
the global and F = F = '^^Tie{n) in the local case, with current conservation resulting from 
TrT; = 0. Conserved vector currents are related to F = — F = Tie{n) in the nonsinglet 
case and to F = — F = e(n) in the singlet case. 



6. Derivation of normal operator 

The Wilson-Dirac operator X/a (with hermitean 7-matrices in 4-dimensional euclidean 
space and < r < 1) is given by 

^ = S E VlV, + m + V 7,(V^ - Vt ) (6.1) 

where (V^)„'„ = 6n'n - Ui,n^n',n+[i (which implies Vj^V^, = V^Vj, = + Vj,). For the 
operator X one has 75-hermiticity, 

Xt = 75X75 , (6.2) 

however, in the presence of a gauge field (with [V^, V^] 7^ and [Vj^, V^] 7^ for 7^ 
and thus [X\ X] 7^ 0) X is not normal. 

To derive a normal and 75-hermitean operator D one needs to start from 

H = 75X (6.3) 

which, being hermitean, in contrast to X has well defined spectral properties. The strategy 
then is, instead of X = 75-/!, to consider 

D = ^^E{H) + C (6.4) 

with some general function E{H) and some constant C, and to determine those quantities 
by imposing the necessary conditions. Requiring 75-hermiticity (|2.2| ) of D it follows that 
E{H) must be hermitean and that C must be real. Since with 

H(j)i = ai4>i (6.5) 

(where ai is real and the (pi form a complete orthonormal set) one has the representation 
E{H) = J2i E{ai)(j)i(j)] , hermiticity of E{H) simply means that E{a) considered as a 
function of a real parameter a must be a real function. 
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From the requirement of normality of D we obtain the condition 

Yi^,E{Hf] = Q . (6.6) 

To satisfy this condition is the central point. Wishing to get a general solution, one must 
require E[H)'^ to be independent of H. Though being inevitable, this is clearly quite 
drastic. It means that E{HY should be a multiple of the identity 

E{Hf = pH , (6.7) 

or that E{a) = ±p with some constant p which, without restricting generality, we can 
take to be positive. In order to keep the properties of E{H) as close as possible to those 
of H we further require E{a) to be nondecreasing and odd. This fixes the signs and we 
end up with 

E{a)=pe{a) (6.8) 

where e{a) = ±1 for a^O. If all 7^ this is already the solution and E{H)/p is just 
the function H/ V H"^ of Neuberger If a; = occur we have to specify e(0). It would 
be tempting to take the value zero for this,0 however, because of ( |6.7| ), i.e. the necessity 
to keep the procedure independent of iJ, this is definitely not possible. Thus one has 
to choose either +1 or —1 for e(0), and one must decide for one of them since no H- 
independent criterion for selection appears available. The oddness of the function, which 
then is violated at a = 0, can be recovered by doing independent calculations for each of 
the two choices and taking the mean of the final results. In Section 7 we will see that in 
terms of counting eigenvalue flows this procedure has a natural equivalent. 

To fix the constant C in ( |6^ ) we note that, because p~^'~^^E{H) is unitary, the spectrum 
of '-)^E{H) is on the circle with radius p and center at zero. Thus to get the appropriate 
spectrum of D we put C = p and get 

D = p{l+^,e{H)) . (6.9) 

It appears important to emphasize that, by the necessity to satisfy ( |6.7| ), one cannot 
escape simultaneously arriving at the Ginsparg- Wilson constraint ( p.5|) on the spectrum 
and at the Neuberger form ( |6.9| ) of the Dirac operator. In addition ( |6.7| ) unavoidably 
produces the somewhat dehcate situation with the choice of e(0). 

To complete the derivation the occurring parameters are to be fixed. The continuum 
limit in the case U = 1 with the representation (V^)pp = 1 — e"*^'''^ indicates that 
masslessness requires m < and one gets p = \m\/a. Further, it is also known that 
to avoid effects of doublers on the finite lattice one needs m > — 2r. A choice with major 
analytical simplifications is —m = r = 1. It should be noted that X enters ( p.9| ) only 

^ Which, by the way, would lead to T{u, v) = ((u — p)^ + w^) ((u — p)^ + v'^ — p^) = instead of ( |3.5D . 
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up to a positive constant factor, so that, for example, using X/a instead of X would not 
change anything. 

Since real eigenvalues occur only at and at 2p, the sum rule ( |2.10| ) for the opera- 
tor (13) reduces to A^+(0) - A^_(0) + N+{2p) - N_{2p) = 0. Similarly ( p7[3|) becomes 
Tr(75L)) = 2p{N+{2p)-N_{2p)). Combining these two relations one has A^_(0)-A^+(0) = 
(2p)~^Tr(75Z}) and inserting the particular form (|6.9| ) of D into this one gets 

iV_(0)-iV+(0) = iTr(e(if)) (6.10) 

for its index. Also the eigenvectors of D and of H can now be related in detail. For 
this we note that ( ^.3| ) in terms of ■y^e^H) becomes '~^^e{H)fk = {^k/p — l)/fc, so that 
parametrizing complex eigenvalues as = p{l + e*"^*) with < < vr and remembering 
( p.8|) we get the eigenequations of e{H) 

e{H)if = TCfc/f for = 




6(i/)/f) = ±/W with /f) = _(e-^'=/Vl'^±e^^'=/Vf) • (6.11) 
These equations are to be compared with 

e(iJ)0i = e(aO0i (6-12) 



which results from ( [6.51 ). Obviously the vectors in ( |6.11D are linear combinations 



i± 



^i^V/^'' where = (pi for e{ai) = ±1. By the properties of D derived here and the 
general structure obtained in Section 2 it is guaranteed that task to find the coefficients 
b^j^^ has a solution. 



7. Relations for spectral flows 

The studies of the fiows of eigenvalues of H with m can be justified on the basis of 
( |6.10| ) which in terms of N^^ and , the numbers of positive and negative eigenvalues 
of -ff, in the absence of eigenvalues zero of H reads 

iV_(0)-iV+(0) = ^«-iV^) . (7.1) 

The crossing of zero of an eigenvalue which occurs at some m is connected to a change of 
the difference of the numbers of positive and negative eigenvalues by +2 or —2, respec- 
tively, depending on the direction of the crossing. Therefore the net number of crossings 
is related to the index of D. 
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To include also zero eigenvalues of H in these considerations, one has to note that in 
the very moment of crossing a positive (negative) eigenvalue has disappeared, however, a 
negative (positive) one has not yet appeared. Still using ( [7. 1|) then agrees with the notion 
that the index in that moment has only changed by |. With this understanding ( |7.1| ) 



is no longer equivalent to ( 6.10 ) in which e(0) = is forbidden by ( |6.7|) . However, the 
analogue of the procedure described in Section 6 (of working with the mean of the choices 
e(0) = +1 and e(0) = —1), in the case of counting flows is seen to lead to the same result 
as the counting at the crossing point mentioned above. Thus the latter appears valid and 
natural. 

To investigate properties of the flows of eigenvalues of H analytically we first derive 
some relations. Multiplying ( |6.5| ) by (pj'-f^ one gets (j)\^^H(j)i = ai(j)\^^(j)i and summing 
this and its hermitian conjugate one has 0; {75, H}(j)i = 2ai(j)\'~)^(j)i. From this by inserting 
( |6.3| ) with ( |0| ) one obtains 



ai (j)l'j5(l)i = gi{m) + m with giim) = (7.2) 



(where < gi{m) < 8r since ||V^0/|| < ||(V^ — l)0z|| + = 2). Next, abbreviating 
(dai)/ (dm) by ai, we note that 

d {(j)lH(j)i) X . ; t H-^ , H-i ^ , d {(j)](f)i) 

— = H(j)i + H(f)i + (f) H(j)i = (f)l -f5(f)i + ai — (7.3) 

dm dm 

which means that we have 

ttz = 0^750; • (7.4) 
Combining (|7^) and ( [7.4|) we get the differential equation 

aiai = gi{m) + m (7.5) 

which can be readily integrated to give 

rm rm 

af {m) = of {rrib) + 2 / dm'(m' + (^^(m')) = a;^(mf,) + m^ — m^ + 2 / dm! gi{m') (7.6) 

in which particular solutions are determined by the choice of af{mf,). 

To get an overview of the set of solutions we note that because of H ^ m75 for m — » ±00 
one gets gi{m) for m —>■ ±00. Therefore, since gi{m) is nonnegative, the equation 
m + gi{m) = has at least one solution mo < 0. If this is the only one we choose nib = ^o- 
Because then /^^ dm'{m' + gi{m')) > for all m it becomes obvious that we can freely 
choose af{mQ) > 0. In this way we get all solutions allowed by ( |7.4| ), which requires di 
to be finite (the forbidden solutions can be conveniently seen by choosing af{mi,) = 
and nib 7^ "^o)- The extension to the general case, where one has to deal with 22; + 1 
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solutions of m + g{m) = with m2z < m2z-i < • • • < '"^i < "'^o, is straightforward. Then 
among the rriy with even y one has to equate that to which leads to the lowest value 
of J^" dm'(m' + giijn') for some fixed m (or one of those in case of degeneracy). 

We thus have a complete specification of the solutions. Clearly all solutions obtained 
show the asymptotic behaviors af(m) — >■ m? for th? oo. It is seen that the points 
iriy with even y determine the characteristic features. If of {my) > there is a minimum 
of the solution +\Jaf{m) and a maximum of the solution —^af(m) at that point. If 
af{my) = then +^Jaf{m) coming from above continues as —\J af{m) below the zero, 

and analogously —^af{m) from above as +^af{m) below, i.e. one gets two solutions 
which cross zero at that point. For the square of the derivative at the crossing point 
(using af = {cxiaiY /of and ( |7.5| )) one obtains 

d;f(m) — > 1 + gi{m) for m — > iriy and af{my) = (7.7) 

which shows that gi{m) in general will cause deviations from the chiral value 1. The 
solutions of the differential equations describe the possibilities for flows which occur. 
Which values of af{mb) and which signs of ±^Jaf{m) are selected and what the detailed 
properties of the function gi{m) are depends on the eigenequation ( |6.5| ) which in the 
present context no longer appears directly. 
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